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Abstract 

The s-channel annihilation of proton and antiproton into a neutral pion and a real 
or virtual photon followed by lepton pair emission is studied. Such mechanism is 
expected to play a role at moderate values of the total energy ^/s, when the pion is 
emitted around 90° in the center of mass. A fair comparison with the existing data is 
obtained taking into account scattering and annihilation channels. The cross section 
is calculated and numerical results are given in the kinematical range accessible in 
the PANDA experiment at FAIR. 



1 Introduction 



The antiproton-proton annihilation in various leptonic and hadronic channels 
will be intensively studied by the PANDA collaboration [1] at the complex 
accelerator FAIR (Darmstadt), where it is planned to accelerate antiproton 
beams with momentum from 1.5 to 15 GeV/c and average luminosity C = 
1.6 ■ 10^^ cm~^ s~^ [2]. The pp annihilation into leptonic channels is very 
challenging, as the hadronic background is higher by orders of magnitude. In 
Refs. [3,4] the annihilation reaction p + p — 7* — )■ e"'' + e~ was investigated 
with the aim to extract of electromagnetic form factors of the proton in the 
time-like region. In this paper we focus on the annihilation reactions 
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p + p ^ 7 + 7r°, 



(1) 
(2) 



which contain also interesting information, in different aspects. Through the 
emission of one virtual photon between the proton-antiproton and the lep- 
ton pair, reaction (1) constitutes a unique opportunity to determine time-like 
proton form factors in the 'unphysical' kinematical region, i.e., under the 
threshold of p + p annihilation into a lepton pair (i.e. for < < 4Mp, where 
q"^ is the momentum transferred to the proton) [5]. Related processes are the 
pion electroproduction on the proton e~ + N ^ e~ + N + tt, which is under 
intensive experimental and theoretical study at MAMI and JLab [6] and the 
pion induced lepton pair production tt + N ^ N + i~ first investigated in 
[7] . Pion scattering on nucleon and nuclei is part of the experimental program 
of the HADES collaboration at GSI [8]. The reaction (1) was studied in Ref. 
[9] near threshold. In Ref. [5] the calculation was extended introducing a larger 
set of diagrams and studying the sensitivity to different parametrizations of 
proton FFs. It was also suggested to extract information of the time-like nu- 
cleon axial FF, through the production of a charged pion from the reaction 
p + n — )■ e+ + -h 7r~ . As no data exist for the nucleon axial FF in the time- like 
region, predictions have been built on analytical extension of models which 
reproduce the data in the space-like region [10]. 

The reaction (1) has also been investigated in backward kinematics, in the re- 
gion of large momentum transfer, with the aim to get information of transition 
distribution amplitudes [11]. In Ref. [12] the process pp ■jir^ was considered 
in the frame of the Generalized Distributions Amplitudes approach. Neither 
the Regge behavior of the scattering amplitude nor the resonant character or 
annihilation amplitude was contained in those papers. 

The reaction mechanism which allows to access the electromagnetic struc- 
ture of the proton is illustrated in Figs. la,b and will be denoted below as 
'scattering' exchange channel. In forward (backward) kinematics the Regge 
character of exchanged nucleon reveals itself. In these kinematical conditions, 
the exchanged nucleon is close to mass shell. Therefore the vertex 7*^*-' * pp 
can be described in terms of two phenomenological quantities Fi 2 which 
can be interpreted as proton form factors (FFs) which describe the partonic 
structure of the proton. This statement is valid in the kinematical region of 
near forward (|t| ^ s) or near backward (|-u| ^ s) scattering. In the region 
\t\ ~ |m| ~ s 3> Mp these amplitudes are (unknown) functions of both kine- 
matical variables. They have the form [M"^ / sy' (j){t / s) , where the exponent n 
is determined by quark-counting rules and are not related with the Dirac and 
Pauli form factors of the proton. A similar behavior is expected for the vertex 
of pion interaction with nucleons. 

In coUinear kinematics (near forward scattering and near backward scatter- 
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ing) the amplitude of the scattering channel dominates, as the contribution 
of the annihilation channel is suppressed by the phase volume factor \t\/s or 
\u\/s. At emission angles near 90°, (center of mass reference frame implied) 
the cross section is dominated by the amphtude of the the 'annihilation' vector 
meson exchange mechanism (see Fig. Ic). Intermediate states such as nucleo- 
nium {pp bound states), vector and scalar mesons (including radially excited 
meson states) can in principle contribute. In Rcf. [13] the creation of narrow 
resonances was discussed. We focus our interest in the energy region outside 
the resonance production. Heavy vector meson intermediate states such as 
uj (1450) and uj (1650) mesons play an important role in limited kinematical 
range, where the total energy is close to their mass, due to Breit-Wigner char- 
acter of relevant amplitudes the states with (1450) and cj(1650). Outside 
this region, they are suppressed by form factors since they are more extended 
objects. 

In this paper, we focus on the processes (1), (2), in the kinematical region 
where the s-channel is expected to dominate (large angle emission of the pion 
and moderate values of the total energy squared, s). The characterization of 
this mechanism (Fig. Ic) is important not only to disentangle the information 
on proton form factors (which are accessible through the 'scattering' mecha- 
nism) Fig. la,b) but also to gain information on the properties of the vector 
mesons. Here the parton (quark) structure of proton turns out to be essential. 
We will limit our considerations to cj-exchange in the intermediate state. It 
is known from previous literature [14,15] that the largest anomalous vertex is 
pujTT, as it has the largest quark coupling. 

The interaction of the vector w-meson with the nucleus in the vertex upp, 
which contains information on the strong proton and meson couplings, can 
be described in terms of proton vector form factors, in the time-like region 
of momentum transfer squared. They have, in principle, complex nature. The 
vertex u nj* can be described either by a phenomenological parametriza- 
tion, or through a triangle vertex. Such triangle vertex, calculable in frame of 
Nambu-Jona-Lasinio model, is relevant also to the transition pion FF [17-19]. 
Recently, large interest has been raised in this field by the new data from 
BaBar on pion FF at large momentum transfer [20]. 

The calculation will be compared to the unpolarized cross section and to the 
angular distribution for p + pj + n annihilation, which was measured in the 
region 2.911 GeV < ^/s < 3.686 GeV by the Fermilab E760 Collaboration 
[21]. 

The plan of the present paper is as follows. In section 2 we calculate the p + 
p — )■ 7 -|- 7r° process for the annihilation channel: the kinematical variables are 
defined, the phase space and the matrix element of the reaction are calculated. 
In section 3 the case of virtual photon is considered and double differential 
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Fig. 1. Feynman diagram for p + p — t- e+ + e + 7r° (a) and (b) in f-channel ('scat- 
tering'); (c) in s-channel ('annihilation'). 

cross section is calculated. In section 4 a realistic parametrization for the 
vertices and the coupling constants are presented. In section 5 we compare the 
calculation to the data in condition of real photon emission around 90°, and 
predictions are given for the double differential cross section for the reaction 
p + p 7* + tt'^. In section 6 (Discussion and Conclusions) the domain of 
applicability of the present model is discussed, as well as its interference with 
the scattering channel exchange mechanism. 



2 Formalism for p + p — > 7 + tt*^ 

Let us consider first the process of proton- antiproton annihilation into a real 
photon of momentum k (with k^ — 0) and a neutral pion: 



Annihilation and scattering exchange channels play both a role, which will be 
discussed below. First the kinematics, which is common to both mechanisms, 
is discussed. 



2. 1 Kinematics 

The calculation is performed in the center of mass system (c.m.s.) (Fig. 2). 

Let us fix the momenta parametrization in c.m.s (in the formp = {pq,P) where 
p is the 3- vector) as: 



p (p_) + p (p+) ^ 7 (A;) + TT (p^) . 



(3) 
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Fig. 2. Schematic view of the reaction p+p ^ 7* + 7r in the center of mass system. 

where E{/3) is the energy (velocity) of the initial proton of mass Mp, ko is the 
energy of the emitted photon and E^^ is the pion energy. Then we can calculate 
the necessary scalar products: 



2{p^q)=AE' = s, 2{qk)=s-Ml 2 = (1 ± /Jc.) , (4) 

where and = cos {j-k,p^ is the cosine of angle between the momenta of the 
initial proton and the produced pion. 

For convenience let us introduce the invariant variables: 



s = {p++ P-f = {k + p^f = g^ u = {p+ - p^f = {p_ - kf , 
t^{p. -p^f ^{p+-k)\ s + t + u^ 2Ml + Ml 

then 



2{p+p_) = s-2Ml 
2{p+p^) = -u + Ml + Ml, 
2{p+k) = -t + Ml 

The phase volume of this reaction 



2{p.p^) = -t + M^ + Ml 
2{kp^) = s-Ml 
2 {p_k) = -u + Ml. 

has the standard form: 



= (27r)' — ^ ^^(5^ (p +p_-k- (5) 

^ ^ {2nf 2ko {27rf 2E^ ^ ^ ^ 

2H^s 2HH ^ 

where the final particles mass-shell (5-functions restricts the energies and the 
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Fig. 3. Feynman diagram for the annihilation reaction p + p — )• 7 + tt^, through 
vector meson exchange. 



moduh of the momenta to: 



2 



k 



and is the pion velocity. The cross section can be calculated using the 
following expression: 

dcT^^jE \M\'d<i>„ (7) 



spins 



where 7 = (1/2) (s - 4M^) = s/3/2 is the invariant flux. 
2.2 p + p^7 + 7r° annihilation through ui meson 



Let us consider the annihilation channel through vector meson exchange of 
the reaction 

P(P-) + P{P+) ^ V{q) ^ 7(/c) + 7r°(pO (8) 



as illustrated in Fig. 3. 

The matrix element for the exchange of intermediate vector meson V can be 
written as: 



where V = cu,... My and Fy (g^) are the mass and the total width of 
intermediate vector meson V with momentum q {(f — s), e'^ is the polarization 
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vector of emitted real photon (i.e., (ek) = 0). In principle the total width of the 
intermediate vector meson depends on g^. This dependence, which is function 
of the produced particle momentum, is particularly important in case of wide 
resonances. It is introduced to suppress unphysical contributions outside the 
resonance region [16]. In the present case, as we are far from the resonance 
mass and we consider u meson which has a small width, we neglect the 
dependence of the width (i.e. s ^ M^, Mt^Ft^). 

The quantity Gv-n-y {(f-, k"^) describes the F (g) — >■ tt (ptt) 7 {k) vertex and it is 
given in section 4. The current related to the ppV vertex J^'> (q^) has the 
form: 

where the vertex F^^^ is parameterized through two (strong) form factors Fi^^ 
in the following way: 

- (e) 7« + ^i^f ^ (f) (11) 

where cr^'^ = (i/2) (7'*7'' — 7*^7'^), and the explicit expression for F^^ is given 
in section 4. 

The region of our interest is far outside the region of narrow resonances which 
is defined as \^/s — Mr\ ~ Fj? <^ ^/s, where Mr(Fr) is the mass(width) of 
the resonance. In the region close to the resonance production the effects of 
formation of bound states is very important and can drastically change the 
value of the cross section [13]. This effect appears in a very narrow range and 
it is not necessary to take it into account in our case. 

In principle, all vector mesons may contribute in the intermediate state. Let us 
give the reasons why the contribution of an intermediate w-meson is dominant 
for this process, with respect to p-exchange. The value of the coupling Vpp is 
9tpp/ (^tt) = 20 » g'^pppl {An) = 0.55. For the Vn^ vertex, the couphng of the 
vector meson with the quark is relevant. The coupling constant of the u and p 
mesons with the valence quarks in the proton or pion are quite different. The 
w-meson, being isotopic singlet constructively interacts with the quarks, unlike 
the p-meson which is isovector and its interaction with quarks is destructive. 
That is the reason for g^^,^^ ~ ^Qpqq- This is a large factor, as the square of this 
constant enters in the amplitude. The annihilation channel contains the vertex 
of interaction of two vector mesons with the pion. This vertex has anomalous 
nature. Therefore, vector mesons should be of different species (for example if 
p meson is in intermediate state, then cu-meson appears in the final state). As 
the decay width of p-meson is about ten times larger that the one for a;-meson 
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we conclude that the main mechanism of the process under consideration is: 

P + p^UJ^TT + p (12) 



with the subsequent transition of the p-meson into a lepton pair or a real 
photon [23]. 

About p', uj' and other radial excitations of p and cj, an additional suppression 
is expected in the considered kinematical region, as form factors decrease 
faster, describing more extended objects. 

Therefore, wc consider only the contribution of the u meson. This simplifies 
also the analytical calculation, as <^ F"^ and can be neglected (see Ref. 
[24]). The vector current takes the simple form: 

J^^^(e)^F'^v{p_)^^u{p^). (13) 
Using the relation 

e^.apce.a^sq'k'^q'k' J2 J'" J*" = 25^^ (2 - /^^ sin^ 0,) \F^f , (14) 

spin 

we obtain for the differential cross section 



, \ ann 

da \ 



where the — s dependence is contained in uo, which has the dimension of a 
cross section. The total cross section as function of is: 

(Tanniq^) = 87rao(?') ( 1 - ^ ) . (16) 



The approximation of w-meson in the intermediate state is well beyond the 
precision on the final cross section, that is evaluated to be ^ 10%. This estima- 
tion takes into account the accuracy of SU (3) symmetry and the contributions 
of the neglected mesons. 
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p(pj yoo 
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P(P ) TZ%p ) 

(a) 



p(pj 5T°(P^) 



► •wv\/v\/\/v- 

P(P,) Y(k) 



(b) 



Fig. 4. Feynman diagram for the reaction p + p ^ j + tt , (a) in t and (b) in u 
channels. 

2.3 t and u channel exchange for p + p ^ 7 + 7r° 



Let us consider the diagrams illustrated in Fig. 4. The corresponding matrix 
element can be written as: 



lb-_ 



KM 



T„ —'jQ 



u 



M2 



u ip+) e", 



(17) 



where e is the proton charge, g^NN is the pseudoscalar coupling constant of the 
pion with the nucleon (the pseudovector coupling has momentum suppression 
[25]), qi = p+ — k, q2 = p+ — Pn and the proton electromagnetic vertex = 
(/c^ = 0) is parameterized in terms of form factors Fi^2' 



r, {k') = F, {k') 7, + {k') a,^k\ (18) 

which, in the real photon limit, coincide with the static values: i^i=l, ^2=1. 798. 
The differential cross section is: 

^<^- = ^EI-^-l'^r2 (19) 



where the phase space can be written as: 



d^kd 



cir2 = (27r)-^- ^5{q'-m'^) = —kodkod cos eS[{p+ + p^ - k)' - M'] 

^dcos9. (20) 



2 

m 



IGtts 

We used the relation (p_,_ + p_ — k^ — — s — 2^/skQ — ml. 
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The differential cross section for t and u channel, and their interference be- 
comes: 

1 <j — 777^ 

da = -^^J2\Msc\'dcose, (21) 

with 



= (t) + 2Fi (t) F2 (t) + (t) , 
F„ = (^) + 2F^ (u) (u) + ^F^ (u) . 

The effects of strong interaction in the initial state interaction which comes 
from exchange by vector and (pseudo) scalar mesons between proton-antiproton 
are really essential here, in the scattering channel. They effectively lead to the 
Regge form of the amphtude of the scattering channel. Thus the t and u dia- 
grams are suppressed by adding a general Regge factor Rt, Ru, in the following 
form: 

«(*)=(-) = (23) 



s 



fl(«)=(-J = - + (24) 



where rag/n ~ 0.7 and sq ci:; 1 GeV^ [26] can be considered fitting parameters. 

This Regge form of amplitude incorporates in principle infinite number of 
resonances, (i.e. A (1232) and others). As for excited resonances like N* (1440) 
they belong to daughter Regge trajectory which contribute power suppressed 
value compared to the leading one. This is also included in the estimation of 
the claimed 10% error. 



3 Formalism for the reaction p + p ^ + e + tt^ 
Let us consider the annihilation channel of the reaction: 

P (P-) +p{p+)^Y (k) + vr° (p.) ^ e+ {k+) + e- (k^) + 7r° (p^ (25) 
where the four momenta of the particles are indicated in parenthesis. 
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We focus here on the mechanism of annihilation through uo meson exchange. 
The general expression of the matrix element, corresponding to the diagram 
of Fig. Ic, for the exchange of a vector meson V , is: 



where e = -\/47ra is the elementary electric charge (a fa 1/137). The lepton 
electromagnetic current has a standard form 



j:^u{k^)Yv{K), (27) 
and the electromagnetic current related to the ppV vertex is 

j; = ^;(p_)rxp+), ri5 = Fr(g^)7^ (28) 



The currents obey gauge invariance: kyj^ = q^iJ^ = 0. Again, only a; ex- 
change is considered and F2 is neglected. 

The cross section of the process (25) can be written in the standard form 



dc7 = j-Jj2\Mfd^s, (29) 

spin 

where is the phase volume of the process: 



^ ^ ^ ^ (27r)'2^+(27r)'2E_(27r)'2E^ 

(27r) ^ 

^^54 _ ^ (30) 

where E± and £^7r are the energies of the electron, the positron and the pion 
in the final state. 

Since our aim is to calculate the angular distribution and energy spectrum of 
the final pion, we insert the unit integration 



J d^k 5^ {k-k+- A;_) = 1. (31) 
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Performing the integration over the final leptons momenta: 



spin \ J 



(32) 



sfnn 

where 



one obtains the following expression for the double differential cross section: 











■(^2 _ + M2r2 





= " " 48,^2fc2 (^') (34) 

A;^ = s + - 2£;^ > 4m^, 

where we put A;^ = in the vertex um^* since we consider the case of small 
virtuality. The energy of the final pion distribution can be obtained by the 
integration over the angle: 



dE. = J ''^dE^d^. = ^°^^^T2?^^ ) V - 



f5^E, = jEl-Ml. (35) 



4 Models of Form Factors and Vertices 



In principle, one should include initial state interaction (ISI). ISI strongly 
modifies the cross section, at near threshold energies. It is evaluated introduc- 
ing a multiplicative factor C, the Coulomb factor, which takes into account 
multiphoton exchange: 

C = — — r, X = — , 36) 
— 1 V 



where v is the incident relative velocity in c units. In the high energy range 
considered here, this factor is close to unity. In our case, this factor increases 
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the cross section, as the interaction occurs between opposite charges. For ex- 
ample, at ^/s = 2 GeV, v = 0.34 c, C = 1.6. As far as meson exchanges are 
concerned, let us remind that in pp annihilation, one pion exchange is sup- 
pressed by the factor Mp^/s which becomes more and more important at high 
energies,the ou meson is included in the ouNN form factor and heavier vector 
meson are suppressed for the reasons given above. 

Following [27], the form factor F{^2 (?^) can be parametrized as: 



with normalization F^{M^) = 1. A^^ = 1.25 GeV is an empirical cut-off. We 
took = 782.65 ± 0.19 MeV. 

This parametrization has been built for the space-like region of momentum 
transfer. In the time-like region, one should consider an extension of this for- 
mula which obeys analyticity, in particular introducing an imaginary part. 
This procedure has been suggested in literature for electromagnetic form fac- 
tors of the miclcon in Rcf. [28], for p- meson [29] and for 04 meson in Rcf. 
[30]. Such paramctrizations are very reliable when data exist to constrain the 
parameters. This is not the present case, and such procedure would add more 
uncertainties. 

The quantity Gi^t^^* is related to the matrix element of the transition U!{q) 
TT{p)j{k) (Fig. 5) and it is in principle momentum dependent: 



M{uj ^ 7r7) = (q\ k') e^^^pq^k^ey^, (38) 

where we introduced the polarization vectors of the photon and of the vector 
meson. 




Fig. 5. Triangle diagram for V tt^. 
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It can be calculated in the following [22] : 



2 a M 



7r3 



(39) 



where we used the Goldbcrgcr-Treiman relation g-Kuu/iT^u = 1/-^?^ — 
93 MeV [23] is the pion decay constant. The constants gumu and gTmu corre- 
spond to the coupling of uj and tt mesons with the light li-quarks in the loop. 
We take the following values [31]: g^juu — 5.94 and g-j^uu — 2.9. 

/ {q^, k"^) is the internal quark loop integral which will be derived in the Ap- 
pendix A. Due to the requirement of absence of real quark in the intermediate 
state, the quantity / (g^, k"^) must be real: 



ml 



In' 



-In' 



mi 



mt 



, e > ml, q^ > m%AO) 



where mu is the constituent quark mass. We verified that our results for real 
and virtual photons agree in the limit k^ m^: 



I (s, k' 



In order to check the consistency of the values of the parameters, we evaluate 
the width of the radiative decay ou — >■ ir^: 



r(a;^A) = — 1-^ « 550 keV. (42) 



192 F2 TT^ V M2 

Thus the decay branching is equal to: 



BR (u ^ 7r%) = ^ 6.5%, (43) 

which is in a fair agreement with the value: BR^^^' {u — )■ 7r°7) = (8.28 ± 0.28) % 
[23]. The values of the constants used in the calculation arc = 8.49 ± 0.08 

MeV, Br{uj n-f) = (8.28 ± 0.28) • lO'^, r^_,^^ = 0.70297 MeV, = 20. 

Alternatively, one can choose a phenomenological parametrization for the ver- 
tex u! — >■ Ti'y^*\ based on monopole dependencies on q^ and k^: 
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where, the constant Gj^tti (0, 0) is derived from the radiative decay a; — >■ 717 
[32]: 



r (a; ^ 717) 



24 




(45) 



5 Results 

No experimental data exist on the reaction p + p e+ + + n^, but the 
reaction p + p — )■ 7 + 7r° was measured in the region 2.911 GeV < y/s < 
3.686 GeV by the Fermilab E760 Collaboration [21] and data exist on the 
cross section and the angular distribution. It is in principle possible to adapt 
the present model to 7 production, and compare the calculation to the data 
in the angular region around 90°, where the present mechanism is expected to 
be dominant (Fig 6). 

A good description with the existing data can be obtained, for parametrization 
(39) of the coupling constant gv^-y, taking into account an additional coupling 
factor of 



due to the interaction of the vector meson with the constituent quarks in 
the vertex Vuu. A good agreement with the data is also obtained using the 
phenomeno logical parametrization (44) and A^^ = 1.25 GeV. This is shown in 
Fig. 7, where the results on the cross section integrated in the range | cos^jrl < 
0.2 are reported and compared with the present calculation in s-channel. 

In Fig. 6 the angular distributions are shown for four values of the total c.m.s. 
energy: (a) = 2.975 GeV, (b) = 2.985 GeV, (c) = 3.591 GeV, 
(d) y/s = 3.686 GeV, together with the results of the calculation. A good 
agreement for the s-channel calculation, Eq. (15), is obtained in the angular 
region around cos^tt ~ 0, because, as expected, this mechanism is applicable 
at these angles. 

In order to reproduce the data in the full range where they are available, we 
have to consider also t and u exchanges. Here we focus to s-channel properties 



A^, 



nl/3 



K + {<f - Kf 



(46) 
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Fig. 6. Angular distributions for different values of the c.m.s. energy: a) -^/s = 2.975 
GcV, (b) ^ = 2.985 GeV, (c) ^ = 3.591 GcV, (d) ^ = 3.686 GeV. The data 
are from [21], the line is the result of the model (see text). 

and we apply a simplified procedure, where we apply the s, t and u channel 
in separated kinematical regions, hence neglecting their interference. 

From the data it appears that the angular distribution flattens in a range 
around cos^ ~ 0, which depends on s. In this central region we apply the 
s— channel result, Eq. (15). In the forward and backward regions, we decon- 
volve in Eq. (23) the u and t contributions and apply the suppression due to 
Regge factors. The results are shown in 6. 

For the s-channel annihilation, which is the object of the present paper, the 
differential and total cross sections have been numerically calculated for the 
reaction p + p — )■ e"*" + e~ + tt" and the results are illustrated for — 7 GeV^ 
and for the parametrization (39) of the (xnrj* vertex. 

The bidimensional plot of the cross section as a function of the c.m.s. angle 
and energy of the pion is shown in 8. The angular and energy dependences are 
fixed by the model. The absolute value of the cross section strongly depends 
on the cutoff parameter taken for the description of the uNN vertex, which 
was taken here as for the real photon data. 

In case of virtual photon, the (7(^777 coupling depends on and k"^. The phe- 
nomenological parametrization can be modified adding a monopole depen- 
dence on k'^. For the parametrization (39), the triangle integral has in general 
more complicated form [33]. Neglecting the pion mass, one can use the sym- 
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q2 [GeV^] 

Fig. 7. (Color online) dependence of the cross section. The data and the calcula- 
tions are integrated for |cos07r| < 0.2 for two choices of the G^^t^^ coupling: from the 
parameterizations (39) (red, dashed line) and from monopole dependence (black, 
solid line). 

metry properties of the triangle loop (see Appendix). 

The bidimensional cross section as a function of the pion variables, at fixed 
total energy s = 7 GeV^ is shown in Fig. 8 and the projection on the pion 
energy is drawn in Fig. 9 for the parameterization (39) of the Gy.,^.^ vertex. 
One can see that the cross section has a smooth angular dependence and it is 
larger when the pion energy is larger. 



6 Discussion and Conclusions 

A hadronic model for the annihilation channel in the processes p + p — )■ vr*^ + 

7 (7*) is presented. This reactions is described by two classes of diagrams. We 
focused here on the annihilation type diagrams. The scattering type diagrams 
with a nucleon in the intermediate state were considered earlier, in Ref. [5] 
and Refs. therein. 

We considered w meson intermediate state as the main contribution because 
it contains the same u and d quarks as the final tt", excluding the p- meson 
due to the small pNN coupling as well as higher resonances. We applied our 
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Fig. 8. (Color online) Bi-dimensional plot of the double differential cross section for 
the process p + p — >• 7* + vr'^, as function of and cos 0,^ at = 7 GeV^. 




0.2 0.4 0.6 0.8 1 1.2 1.4 
E. [GeV] 



Fig. 9. Differential cross section for the process p + p — )• 7* + vr'', as a function of 
at g2 = 7 GeV^. 
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calculation in the energy region outside the resonance production. We com- 
pared our calculation to the existing data on pp — )■ 7r°7, and found a good 
description taking into account also t and u channel, in forward and backward 
angular regions. The formulas obtained in this paper describe the annihila- 
tion channel through a;-meson (which gives an almost isotropic contribution 
in the center of mass system) as well as the near forward or near backwards 
kinematics (Rcggeized nuclcon in t and u channels). Wc reproduce the experi- 
mental evidence that the range of the central angular region, where s-channel 
dominates, increases when the energy increases. 

As it was shown in Ref. [5], the contribution of the t-channel can be pa- 
rameterized in terms of two unknown functions F1.2 (t, u) which depend of 
both variables t and u {s + t + u = 2M^) and rapidly decrease in the region 
s ~ |t| ~ \u\ as Fi^2 {t, u) ~ (M^/s)" <^ 1 where n > 1. Let us note that in the 
limits |i| <S s or <^ s these quantities do coincide with the electromagnetic 
form factors of the proton. In these kinematical regions the Regge-factors have 
to be included into the differential cross section. These factors rather strongly 
suppress the cross section. 

Initial state interaction through meson exchange is effectively taken into ac- 
count in our calculation, in the following way: - in the annihilation channel 
it is in principle included in the form factors which describe the vertex ppu, 
which is in principle complex. We should note that in model which we use 
and which was described in [27] this coupling is real. In general it is possible 
to use parametrizations built for the space-like region and make an analytical 
extension to time-like region. However any extension of the parametrization 
requires extra parameters, which can not be presently constrained, and induce 
unnecessary complications in the model. - in the scattering channel the effects 
of ISI are included in the Regge factor of the amphtudes. 

Let us note that in case of 7* production, with subsequent conversion to 
electron positron pair, in the annihilation channel, one can take into account an 
intermediate production of p meson form the ou by the following replacement: 

^ Gpe+e^Grpuu fA7\ 

¥ - Ml + iMpVp ^ ' 



where Gpuu = Gppp and Gpe+e- is determined from the width of the decay 
p — >■ e"*" -|- e~: 

Gpe.e- - , (48) 



with Tpe+e- = 7 • 10-3 MeV [23]. 
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We estimate an accuracy of our calculation of 10% which is evaluated from 
the precision of Nambu-Jona-Lasinio model and the contributions of other 
vector mesons which we do not consider here. 

The present approach can be generalized to all pseudoscalar mesons, tt, 77, 77/... 
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A Loop integral 



In (39), the vertex V (q) ^ tt (p^^) 7 (k) was expressed in terms of quark loop 
integral: 



d^Pq 



I 



iTT^ (p2 _ ^2 + ,0) ((p, - qf -ml + iQ) [{pq - kf -ml + iQ 

1 



1 l-x 

= I dx I dy- -5 

1 - -^xy - -^xz - ^vz - lO 

mi ^ mi mi 



(A.l) 







where m„ is the mass of w-quark which circulates in the loop. The value of this 
integral in case of real photon (i.e. k? — 0) is well-known (see, for instance, 
Ref. [34]): 



ml' ml) ' 



l{q',0)=l2 
where I2 {a,b) was defined in [34]: 



(A.2) 



l2{a,b) 



a — b 



(A.3) 



20 



arcsm 



2v^ 

log - tTT 



X > 



(A.4) 



and r]± = |(1 ± ^/T^^). 



Applying this formula in the triangle Feynman diagram to the quark loop 
in the intermediate state, one should drop the imaginary term itt (Fig. 5). 
Following confinement properties, real quarks on mass shell are not allowed. 

In the process involving a virtual photon, (i.e., k > 0), this integral has a 
more complicated form. However one can find a simple expression neglecting 
the pion mass. Indeed, the expression (A.l) is symmetric with respect to 5^, 
A;^, and M^, therefore, one can use Eq. (A.4) with 
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